Third minima in thorium and uranium isotopes in the self-consistent theory 
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Background: Well-developed third minima, corresponding to strongly elongated and reflection-asymmetric 
shapes associated with di-molecular configurations, have been predicted in some non-self-consistent models to 
impact fission pathways of thorium and uranium isotopes. These predictions have guided the interpretation of 
resonances seen experimentally. On the other hand, self-consistent calculations consistently predict shallow third 
minima. 

Purpose: We study the isentropic potential energy surfaces of selected even-even thorium and uranium isotopes 
at several excitation energies. In order to understand the driving effects behind the presence of third minima, we 
study the interplay between pairing and shell effects. We also investigate the interpretation of third minima in 
terms of di-molecular (cluster) configurations. 

Methods: We use the finite-temperature superfluid nuclear density functional theory. We consider several 
Skyrme energy density functionals, including traditional functionals such as SkM* and a new functional UN- 
EDF1 recently optimized for fission studies. 

Results: The potential energy surfaces for 228 - 232 Tri and 232 U at several excitation energies are presented. We 
also present isotopic chains to demonstrate the evolution of the depth of the third minimum with neutron number. 

Conclusions: We demonstrate that the depth of the third minimum is sensitive to the excitation energy of 
the nucleus. In particular, the thermal reduction of pairing, and related enhancement of shell effects, at small 
excitation energies help to develop deeper third minima. At large excitation energies, shell effects are washed out 
and third minima disappear altogether. According to our calculations, third minima become more pronounced in 
neutron-deficient isotopes such as 226 . 228 Th and 228 - 230 U. 

PACS numbers: 24.75. +i, 21.60.Jz, 27.90.+b, 24.10.Pa 



I. INTRODUCTION 



The phenomenon of nuclear fission is a large ampli- 
tude collective motion in which the nucleus undergoes a 
series of shape rearrangements before splitting into dis- 
tinct fragments. The observables for a fissioning system, 
such as fission half-life and properties of fission fragments, 
are sensitive to the sequence of nuclear configurations 
through which the nucleus is driven on the way to fission 
[IHH . Local minima in the potential energy surface, often 
representing mctastable configurations, can profoundly 
affect the dynamics and timescale of fission. Of particu- 
lar importance are superdeformed fission isomers [1, Q , 
corresponding to the "second mininum" in actinide nu- 
clei, separating inner and outer saddles. Another impor- 
tant class of states are hyperdeformed "third minima" 
predicted theoretically in the early seventies 0, 0] and 
soon afterwards attributed to the resonance microstruc- 
tures observed in the fission cross sections found in the 
light actinides [|| . Continued experimental studies of the 
actinides (9l-[20j inferred the existence of highly elongated 
minima, and its reflection asymmetric structure has been 
supported by the presence of parity doublets 21 1. 

The appearance of third minima around 232 Th has 
been attributed to large shell effects associated with re- 
flection asymmetric configurations corresponding to di- 



molecular stuctures |22| - |25| . with one fragment resem- 
bling the doubly-magic 132 Sn (2fJ 28]. Pronounced third 
minima have been predicted in theoretical studies of the 
thorium and uranium isotopes, especially those carried 
out with the macroscopic-microscopic (MM) approach 
d 0, [M H3, HI, On the other hand, self-consistent 
studies based on the nuclear density functional theory 
[3ll - [35j . as well as recent MM work [3(|, typically find 
a third minimum that is much shallower than that of 
the earlier MM calculations or the empirical barrier fits 
[3?l |38| . This result is puzzling in light of of the accumu- 
lated experimental evidence (resonances in fission cross 
sections, mass and kinetic energy distributions of fission 
fragments, fits to experimental cross sections, moments 
of inertia, and presence of parity doublets) . 

To clarify the situation, we carry out self-consistent 
calculations for eight even-even Th and U isotopes within 
the superfluid, finite-temperature nuclear DFT theory 
(FT-DFT), investigating how the potential energy sur- 
face and the third minima evolve with excitation energy. 
In particular, we seek to isolate the contributions to the 
nuclear energy that may be responsible for third min- 
ima. We review the FT-DFT model in Sec. HI SectionJH 
presents an analysis of the trends seen in potential energy 
curves and shell correction energies across Th and U iso- 
topic chains, finding that deep third minima appear in 
lighter isotopes. We proceed to analyze two-dimensional, 
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finite-temperature potential energy surfaces of 228 > 232 Th 
and 232 U in Sec. IIVI As the excitation energy increases 
and pairing quenches, we actually find a regime in which 
the third minimum is slightly deepened for 232 Th. Fi- 
nally, the conclusions of our work are given in Sec. [Vj 



II. THE MODEL 

To study the potential energy surfaces (PES) as a 
function of the excitation energy E* , we employ the su- 
perfluid FT-DFT theory [39|, H(| in the implementation 
of Refs. [H, Hi], [42j. To solve the finite-temperature 
Hartree-Fock-Bogoliubov (HFB) equations, we employ 
the symmetry- unrestricted Skyrme DFT solver HFODD 

sua. 

To constrain the total quadrupole moment Q20 (elon- 
gation) and total octupole moment Q30 (reflection asym- 
metry, important at outer saddle and on to scission) we 
employ the augmented Lagrangian method [45[ . We also 
apply a constraint on the triaxial quadrupole moment 
Q22 to force the system to break axial symmetry, sub- 
sequently relaxing this constraint to allow the system to 
follow the minimum energy path in the Q22 direction. 
Consequently, the inner fission barrier heights are low- 
ered due to the axial symmetry breaking. 

The finite-temperature HFB equations are obtained 
from the minimization of the grand canonical potential, 
so that the free energy F = E — TS is formally calcu- 
lated at a fixed temperature T. Since the system is not 
in contact with a heat bath, the fission process is not 
isothermal; but since the large-amplitude collective mo- 
tion is much slower than the single-particle motion, it is 
reasonable to treat fission as an adiabatic process (4r| . 
By assuming fission to be an isentropic process, we ex- 
ploit the correspondence between surfaces of free energy 
at constant temperature and surfaces of internal energy 
at constant entropy [U, |U |H| . That is, we calculate the 
free energy for a fixed temperature as a function of the 
collective coordinates, understanding that relative quan- 
tities such as barrier heights will identically match those 
obtained from a calculation of internal energy at fixed 
entropy. This identity has been verified numerically in 
the self-consistent calculations of Ref. [4l| . 

To map the excitation energy of the nucleus E* to the 
fixed temperature T, we use the correspondence 

E*(T)=E g . s .(T)-E e . s .(T = 0), (1) 

where E S _ S _(T) is the minimum energy of the nucleus at 
temperature T. This corresponds well to the excitation 
energy of a compound nucleus [U, [42j] . 

To study the role of shell effects in producing the third 
minimum, it is useful to employ the shell correction of 
the MM approach. Using the Strutinsky energy theorem 
[49j , the self-consistent energy E can be decomposed as 

E = -^smooth + SE sh , (2) 



where -^smooth is a bulk contribution to the energy that 
varies smoothly with nucleon number and 5E sh is a shell 
correction energy. To extract 8E sh from HFB energy, 
we employ the procedure described in Refs. [5(1 HH with 
the smoothing width parameter 7„ = 1.54, 7 P = 1.66 (in 
units of huiQ = 41 /A 1 / 3 McV) and the curvature correc- 
tion p = 10. 

The nuclear interaction in the particle-hole channel 
has been approximated through the SkM* parametriza- 
tion [52j ] of the Skyrme energy density functional (EDF). 
This traditional EDF achieves realistic surface prop- 
erties in the actinides, allowing a good description of 
the evolution of the energy with deformation [53l [54| . 
At E* = MeV, we also performed calculations with 
the recently developed EDF parametrization UNEDF1 
[55| . In the particle- particle channel, we use the density- 
dependent mixed-pairing interaction (56|. All calcula- 
tions were performed with a quasiparticle cutoff energy 
of E cut = 60 McV. The pairing strengths V r o (r = n,p) 
are chosen to fit the pairing gaps determined from ex- 
perimental odd-even mass differences in 232 Th [H3|. For 
SkM* EDF, the pairing strength are V n0 = -273.5 MeV 
and Vp = -334.0 MeV. 



III. THIRD MINIMA IN Th AND U ISOTOPES 

Self-consistent calculations tend to predict either a 
very shallow or no third minimum for 232 Th and 232 U. 
As seen in Fig. [TJ the potential energy curves for 232 Th, 
obtained with several EDFs, each exhibit a gentle down- 
wards slope beyond the second saddle (Eb) — none of 
these models predicts a large third hump. A shallow 
third minimum appears around Q20 = 150 b in our SkM* 
model and the HFB-14 calculations of Ref. [5§ |. This 
minimum seems to be more pronounced in the relativistic 
DFT calculations of Ref. j3y employing PL-40, NL1, and 
NL-SH functionals. In the Gogny D1S model of Ref. [H, 
the minimum is shifted to Q20 ~ 200 b. (Since they did 
not use the same definition of Q20 as in this work, the hor- 
izontal scale for D1S is not precise.) Our UNEDF1 results 
show a local plateau at Q20 ~ 200 b but the third barrier 
is practically nonexistent. The results shown in Fig. Q] 
and Ref. [31J indicate that self-consistent models predict 
a third minimum, or a softness in the PES of 232 Th in the 
region beyond the outer saddle with Q20 ~ 150 — 200 b. 
This indicates that the shell effects responsible for this 
structure are systematically present in DFT calculations, 
but their strength is strongly model dependent. 

The third minimum of 232 Th has been associated with 
a dimolecular configuration, in which one fragment bears 
a strong resemblance to 132 Sn (2(| [23]. The nuclear den- 
sity profiles of 232 Th, 132 Sn, and 100 Zr, calculated with 
SkM*, are displayed and compared in Fig. [2J When 
making this comparison, we followed the methodology 
of Ref. [54|]. Namely, the configuration of 232 Th corre- 
sponds to the local third minimum at Q20 = 165 b, the 
configuration of 132 Sn corresponds to its spherical ground 
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FIG. 1. (Color online) Potential energy curves for 232 Th ob- 
tained with several EDFs: SkM* and UNEDFl (this work), 
HFB-14 and D1S [H. The empirically inferred values 
of the first and second barrier heights Ea and Eb 1371 ]. as 
well as the measured energy of the fission isomer En [5y| , are 
marked. 




FIG. 2. (Color online) (a) Cross section of the total density 
of 232 Th in the yz-plane calculated in SkM* at the third min- 
imum (Q20 = 165 b), compared to the cross sections of 132 Sn 
and 100 Zr densities, (b) Density profiles for 2:!2 Th and the 
1,!2 Sn and 100 Zr fragments along the z-axis. 



state, and the 100 Zr fragment configuration corresponds 
to its prolate ground state with Q20 = 10 b. The re- 
semblance of the left-hand fragment of 232 Th to 132 Sn is 
clearly seen, although the nascent fragments overlap to 
produce the sizable neck seen in the figure. Similarly, 
our calculations result in a dimolccular configuration for 



232 U predicted to consist of 132 Sn plus a slightly deformed 
100 Mo. 

As discussed in, e.g., Refs. [H, 0, H3] , the high like- 
lihood of obtaining 132 Sn-like fragments in the fission of 
actinides can be attributed to the doubly-magic nature 
of 132 Sn. The recent theoretical studies of the asymmet- 
ric fission around 180 Hg [54], [6(| indicate that the shell 
effects at pre-scission configurations associated with the 
deformed fragment play an equally significant role in the 
determination of fission yields. 

A more comprehensive DFT survey of the actinides 
shown in Fig. [3]Ja,b) reveals that the lighter isotopes 
of thorium and uranium, 22 6> 22 8,230 Th and 228,230 U; are 

expected to have deeper third minima. The hypcrdc- 
formed minimum in 228 Th, for example, appears to be 
a dimolecular 132 Sn-plus- 96 Zr configuration. Why does 
this configuration result in a deep third minimum? We 
seek an answer in the underlying shell effects displayed 
in Fig. EM). 
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FIG. 3. (Color online) The potential energy curves (top) 
and total shell corrections 5E Bh (bottom) predicted for 

226,228,230,232 Th (lgft) 228,230,232,234 y ^ 

The shell corrections for the N — 136, 138 isotopes, 
226,228 Th and 228,230^ indicate a strong shell effect at 
Q20 ~ 150 b. For the TV = 140,142 isotopes ( 23 °' 232 Th 
and 232 > 234 U), the total shell correction tends to stabilize 
more elongated configurations, at Q20 ~ 200 b. This re- 
sult is reminiscent of a spherical-to-deformed shape tran- 
sition around N = 58 in the Zr and Mo isotopes [6l[ - 
associated with the lighter fragments in the dimolccular 
picture of the third minimum. 

To see whether this is indeed the case, in Fig. [4] we plot 
the individual neutron and proton shell corrections for Th 
and U isotopes as a function of Q2o- As expected, proton 
shell corrections weakly depend on the neutron number, 
and they all exhibit a minimum around Q20 = 160 b. For 
N = 142, the neutron shell correction shows two min- 
ima: one around Q20 — 120 b and the second one around 
Q20 = 200 b. While the first minimum weakly depends 



4 




100 200 100 200 

Q 20 (b) 

FIG. 4. (Color online) Neutron (solid lines) and proton 
(dashed lines) shell correction energies as functions of Q20 
for Th (left) and U (right) even-even isotopes calculated in 
SkM*. 



on N, the second one is absent in N = 136, 138 isotones. 
It is tempting, therefore, to associate the large neutron 
shell effect at Q20 ~ 200 b with the prolate-deformed 
N ps 60 fragments, and the large neutron shell effect at 
Q20 ps 120 b with the nearly-spherical TV ps 54 fragments. 
Since the maximum of the proton shell effects appears at 
the minimum of the neutron shell effect, the total shell 
correction is sensitive to both N and Q20- This explains 
the shallow third minima obtained in DFT calculations. 

In summary, we have found that our self-consistent 
SkM* model predicts enhanced third minima for the 
N = 136, 138 isotopes of Th and U. Since shell effects 
are thermally quenched, in the next section we turn to 
study PESs of 22 ^ 2 32Th and 232 U as a function of E*. 



IV. EXCITATION ENERGY DEPENDENCE 

To discuss the excitation-energy dependence of the 
third minimum, Fig. [5] displays the isentropic potential 
energy curves at constant entropy for 232 Th at several 
excitation energies. As excitation energy increases from 
E* = to E* = 48MeV, the second fission barrier is 
gradually reduced while a third barrier changes little. 
This deepens the third-minimum pocket. The apparent 
stabilization of the third minimum at intermediate values 
of E* can be attributed to the interplay between pairing 
and shell effects [46|, |62|. Indeed, as discussed in, e.g., 
Ref. [63j . as the excitation energy increases, pairing cor- 
relations are quenched faster than the shell effects. This 
gives rise to a reentrance of shell effect with E* in the 
third barrier region, so that the third minimum becomes 
more pronounced for moderate excitation energies. Be- 
tween E* rs 21McV and E* ps 48MeV, the second bar- 
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FIG. 5. (Color online) Isentropic potential energy curves for 
232 Th computed in SkM* at several values of excitation en- 
ergy (in MeV). The minimum potential energy at given E* is 
always normalized to zero. 



rier vanishes but the extended plateau around the third 
minimum is still visible. 

To illustrate how the energy plateau around the third 
minimum For 228 < 232 Th and 232 U changes with E*, we 
display isentropic PESs in the (Q20,Q3o) plane in Fig. EI 
By constraining the triaxial moment Q22 , we account for 
the effect of triaxiality on inner barriers. We trace the 
lowest-energy pathway from the ground state to the exit 
point of the barrier. 

For all three isotopes, there is a strong preference for 
an asymmetric fission pathway at low energies. As E* in- 
creases, the barrier to symmetric fission lowers substan- 
tially so that the symmetric fission gradually becomes to 
compete with the asymmetric channel. This describes 
the situation seen experimentally: the mass distribution 
of fission fragments in actinides is strongly asymmetric 
at low energies, and the symmetric mass yield increases 
with E* . For example, the experiment of Ref. 6J| mea- 
sured the mass yield for the photofission of 232 Th, re- 
porting that the ratio of symmetric yield to asymmetric 
yield increases from 2% to 10% for a bremsstrahlung en- 
ergy range (corresponding approximately to our excita- 
tion energy) of 15 — 55McV. 

As seen in Fig. [SJ the third minimum in 228 Th is ac- 
tually rather robust - a 1 MeV pocket still remains at 
E* = 21 MeV, where pairing is completely quenched. 
At higher excitation energies, the third minimum dis- 
appears as the symmetric fission pathway opens. The 
two-dimensional PES shows the same evolution for 232 Th 
as seen in the one-dimensional plots of Fig. [5j Namely, 
as pairing is quenched, the plateau around Q20 = 160 b 
deepens into a well-developed third minimum. Again, 
at higher energies the third minimum disappears as the 
symmetric fission channel opens. For 232 U, however, 
our SkM* calculations do not appear to predict a clear 
third minimum anywhere in the range of E* studied. A 
potential-energy shoulder appears at the lowest energies, 
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FIG. 6. (Color online) Isentropic potential energy surfaces in the (Q20, Q30) plane for 228 . 232r ph and 232 U calculated with SkM* 
at E* = 0, 21, and 47 MeV. The static fission pathways are indicated. A constraint on the triaxial quadrupole moment Q22 has 
been applied to minimize the total energy in the direction of Q22- Consequently, the inner fission barriers are lowered due to 
the breaking of axial symmetry. 



so this may be the source of the resonances observed in 
the experimental data presented in Ref. (l8j . 



V. CONCLUSIONS 

This self-consistent FT-DFT study predicts very shal- 
low third minima, or shoulders, in the potential energy 
surfaces of 232 Th and 232 U. Those minima correspond to 
hyperdeformed and reflection-asymmetric shapes. In the 
lighter isotopes with TV = 136 and 138, the third min- 
ima are better developed. This can be traced back to the 
neutron shell effect that reduces the third outer barrier 
at Q20 ~ 200 b at N = 140 and 142. The shallowness of 
the third minimum in 232 Th and 232 TJ is a robust feature 
of many DFT calculations, including SkM*, UNEDF1, 
D1S, and HFB-14 models. 

Our study fully confirms earlier suggestions that the 
third minimum can be associated with a dimolecular con- 
figuration involving the spherical doubly-magic 132 Sn and 
a lighter Zr or Mo fragment. We show that the neutron 
shell effect that governs the existence of the third min- 
imum, and makes the third minimum more pronounced 



in N = 136 isotopes as compared to N = 142 systems, 
is consistent with the spherical-to-deformed shape tran- 
sition in the Zr and Mo isotopes around N = 58. 

The third minimum is found to be quite sensitive to ex- 
citation energy. Our FT-DFT study predicts that third 
minima become deeper at moderate excitation energies, 
where pairing correlations are quenched and shell ef- 
fects become locally enhanced. At large values of E*, 
the conditions needed for the hyperdeformed metastablc 
states to exist deteriorate as the symmetric fission chan- 
nel opens up. 

While the inference of a hyperdeformed fission isomer 
from experimental data does rely on many assumptions, 
the accumulated experimental evidence is substantial and 
should not be considered lightly. Shallow third minima 
(or shoulders) obtained in self-consistent calculations are 
in fact consistent with the observations. The absence of 
a well-developed local minimum in a static PES, a sole 
focus point of Ref. [36( , does not tell the full story. Of- 
tentimes, observed states can be associated with config- 
urations, which do not correspond to a minimum in PES 
pfjl [6o| but are well separated from other states through 
the presence of specific quantum numbers. Also, the en- 
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crgctics of local minima can be impacted by correlations 
associated with symmetry restoration, such as those dis- 
cussed in Refs. [6a Ir37|. 
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